lopscience 


iopscience.iop.org 


Home  Search  Collections  Journals  About  Contact  us  My  lOPscience 


Influence  of  strain  on  thermal  conductivity  of  silicon  nitride  thin  films 


This  article  has  been  downloaded  from  lOPscience.  Please  scroll  down  to  see  the  full  text  article. 
2012  J.  Micromech.  Microeng.  22  045001 
(http://iopscience.iop.org/0960-1 31 7/22/4/045001 ) 

View  the  table  of  contents  for  this  issue,  or  go  to  the  journal  homepage  for  more 


Download  details: 

IP  Address:  131.84.11.215 

The  article  was  downloaded  on  14/03/2012  at  14:04 


Please  note  that  terms  and  conditions  apply. 


Form  Approved 
OMB  No.  0704-0188 


Report  Documentation  Page 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

MAR  2012 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

4.  TITLE  AND  SUBTITLE 

Influence  Of  Strain  On  Thermal  Conductivity  Of  Silicon  Nitride  Thin 
Films 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory, Thermal  Sciences  and  Materials 

Branch, Wright-Patterson  Air  Force  Base, OH, 45433 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Journal  Of  Micromechanics  And  Microengineering,  2  March  2012 

14.  ABSTRACT 

We  present  a  micro-electro-mechanical  system-based  experimental  technique  to  measure  thermal 
conductivity  of  freestanding  ultra-thin  films  of  amorphous  silicon  nitride  (Si3N4)  as  a  function  of 
mechanical  strain.  Using  a  combination  of  infrared  thermal  micrography  and  multi-physics  simulation,  we 
measured  thermal  conductivity  of  50  nm  thick  silicon  nitride  films  to  observe  it  decrease  from  2.7  W  (m 
K)&#8722;1  at  zero  strain  to  0.34  W  (m  K)&#8722;1  at  about  2.4%  tensile  strain.  We  propose  that  such 
strong  strain?thermal  conductivity  coupling  is  due  to  strain  effects  on  fraction?phonon  interaction  that 
decreases  the  dominant  hopping  mode  conduction  in  the  amorphous  silicon  nitride  specimens. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

10 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


IOP  Publishing 


Journal  of  Micromechanics  and  Microengineering 


J.  Micromech.  Microeng.  22  (2012)  045001  (8pp)  doi:  10. 1088/0960- 13 17/22/4/045001 

Influence  of  strain  on  thermal 
conductivity  of  silicon  nitride  thin  films 

M  T  Alam1,  M  P  Manoharan1,  M  A  Haque1,3,  C  Muratore2 
and  A  Voevodin2 

1  Department  of  Mechanical  and  Nuclear  Engineering,  Penn  State  University,  University  Park,  PA  16802, 

USA 

2  Air  Force  Research  Laboratory,  Thermal  Sciences  and  Materials  Branch, 

Wright-Patterson  Air  Force  Base,  OH  45433,  USA 

E-mail:  mah37  @psu.edu 

Received  20  October  2011,  in  final  form  20  January  2012 

Published  2  March  2012 

Online  at  stacks.iop.org/JMM/22/045001 

Abstract 

We  present  a  micro-electro-mechanical  system-based  experimental  technique  to  measure 
thermal  conductivity  of  freestanding  ultra-thin  films  of  amorphous  silicon  nitride  (S^N^  as  a 
function  of  mechanical  strain.  Using  a  combination  of  infrared  thermal  micrography  and 
multi-physics  simulation,  we  measured  thermal  conductivity  of  50  nm  thick  silicon  nitride 
films  to  observe  it  decrease  from  2.7  W  (m  K)-1  at  zero  strain  to  0.34  W  (m  K)-1  at  about 
2.4%  tensile  strain.  We  propose  that  such  strong  strain-thermal  conductivity  coupling  is  due  to 
strain  effects  on  fraction-phonon  interaction  that  decreases  the  dominant  hopping  mode 
conduction  in  the  amorphous  silicon  nitride  specimens. 

(Some  figures  may  appear  in  colour  only  in  the  online  journal) 


1.  Introduction 

Multi-domain  (electrical,  mechanical  and  thermal,  to  name  a 
few)  coupling  is  an  essential  feature  of  micro-  and  nanoscale 
devices  and  systems  [1] .  The  same  is  true  for  other  applications 
in  micro-electronics  [2,  3],  opto-electronics  [4,  5],  flexible 
electronics  [6,  7]  and  energy  conversion  [8,  9].  These  devices 
are  commonly  fabricated  out  of  thin  film  materials,  which 
experience  mechanical  strain  during  deposition  as  well  as 
device  operation.  Mechanical  reliability  of  thin  films  is 
therefore  as  important  as  their  transport  (electrical  and  thermal) 
properties.  While  electro-mechanical  coupling  is  commonly 
exploited  in  microelectronic  and  other  small-scale  systems,  the 
thermal  domain  must  also  be  considered  in  these  applications, 
whether  due  to  intentional  design  or  unintentional  coupling 
(for  example,  excessive  heat  generation  in  chips  [10]). 
Thermo-mechanical  coupling  in  thin  films  is  typically  studied 
in  terms  of  temperature  effects  on  mechanical  properties 
[11-13]  and  not  from  a  thermal  transport  perspective.  The 
fundamental  question  whether  mechanical  strain  can  influence 
thermal  transport  is  therefore  the  core  theme  for  this  study. 

3  Author  to  whom  any  correspondence  should  be  addressed. 


In  this  paper,  we  investigate  low  pressure  chemical 
vapor  deposited  silicon  nitride  (LPCVD  S^NzO  thin  films. 
Silicon  nitride  is  prevalently  used  for  isolation,  passivation, 
etch  masking  as  well  as  structural  and  optical  layers  for 
various  micro-electronic,  opto-electronic  and  micro-electro- 
mechanical  systems  [14-17].  It  also  experiences  very  large 
residual  strain  during  deposition  and  exerts  similar  strains  to 
other  thin  film  materials  attached  to  it,  which  makes  it  a  suitable 
candidate  for  this  study.  In  fact,  the  large  elastic  strain  has  led 
researchers  to  propose  silicon  nitride  substrates  to  tune  the 
electron  mobility  in  microelectronic  devices  [18,  19].  Further 
motivation  for  this  study  comes  from  the  observation  that 
the  role  of  mechanical  deformation  on  thermal  conductivity 
has  been  computationally  explored  by  only  a  few  researchers 
[20-24]  and  experimental  efforts  have  appeared  in  the 
literature  only  very  recently  [25-27]. 

A  brief  review  of  the  state  of  the  art  in  thermal  conductivity 
of  silicon  nitride  is  given  below.  Mastrangelo  et  al  [28] 
measured  thermal  conductivity  of  microns  thick  silicon  nitride 
films  to  be  around  3.2  W  (m  K)-1.  Zhang  and  Grigoropuolos 
[29]  observed  anomalous  thickness  dependence  and  suggested 
that  microstructural  defects  may  strongly  influence  thermal 
conductivity.  Jain  and  Goodson  [30]  measured  the  in-plane 
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Figure  1.  ( a )  Schematic  diagram  showing  a  single  layer  silicon  nitride  three-strip  actuator  and  a  specimen.  ( b )  Scanning  electron 
micrograph  of  a  specimen  integrated  with  a  five-strip  actuator  and  a  thin  film  heater,  (c)  Strain  is  measured  after  post-experiment  fracture  of 
the  specimen. 


thermal  conductivity  of  1 .5  gm  thick  specimens  to  be  about 
5  W  (m  K)_1.  At  the  nanoscale,  Sultan  et  al  [31]  reported 
thermal  conductivity  of  500  nm  thin  films  as  3-4  W  (m  K)-1  for 
a  temperature  range  of  77-325  K.  For  180-220  nm  thick 
low  stress  nitride,  Zink  and  Heilman  [32]  observed  stronger 
temperature  effects  thermal  conductivity  varying  from  0.07 
to  4  W  (m  K)-1  at  a  temperature  range  of  3-300  K.  The 
cross-plane  thermal  conductivity  measured  by  Lee  and  Cahill 
[33]  and  Shin  et  al  [34]  for  less  than  100  nm  thickness  was 
in  the  range  of  0.4-0.7  W  (m  K)-1,  showing  very  strong 
pseudo-size  effects,  which  were  ascribed  to  the  interfacial 
thermal  resistance.  Stojanovich  et  al  [35]  measured  thermal 
conductivity  of  180  nm  thick  freestanding  specimens  to  be 
about  2.1  W  (m  K)_1.  Bai  et  al  [36]  measured  thermal 
conductivity  of  stoichiometric  silicon  nitride  to  be  about  1.2- 
2.0  W  (m  K)-1  for  thicknesses  of  37-200  nm  thick  films. 
It  is  very  important  to  note  that  none  of  the  above  studies 
measure  thermal  conductivity  as  a  function  of  mechanical 
strain.  Even  though  the  specimens  are  typically  freestanding, 
the  mechanical  boundary  conditions  result  in  about  0.2%- 
0.3%  residual  strain  depending  on  the  silicon- to-nitrogen  ratio. 

2.  Device  design  and  fabrication 

Silicon  nitride  is  a  very  brittle  material  with  unusually 
large  tensile  residual  stress,  which  makes  its  thermo-physical 
characterization  challenging.  Common  approaches  are  to 
synthesize  low  stress,  silicon  rich  (SiNx,  x  =  1-1.1)  nitride 
or  to  study  the  specimen  on  a  rigid  substrate.  Because  the 
interfacial  thermal  resistance  is  challenging  to  measure  and 
plays  a  significant  role  in  nanoscale  experiments,  studies  on 
freestanding  specimens  are  desirable.  Characterization  and 
control  of  the  amount  of  stress  or  strain  in  ultra-thin  specimens 
are  also  challenging.  To  obtain  specimens  with  different  values 
of  strain,  one  needs  to  vary  the  deposition  parameters  from 
wafer  to  wafer,  which  is  cumbersome  and  also  changes  the 
composition  and  micro  structure  of  the  specimen.  We  address 
this  issue  with  the  concept  of  ‘self- actuation’  by  exploiting  the 
tensile  residual  stress  in  silicon  nitride.  Here,  we  design  both 


actuator  and  the  specimen  in  the  same  film  plane  but  increase 
the  cross-sectional  area  of  the  actuator  portion  compared 
to  the  specimen.  This  is  shown  in  figure  1(a) ,  where  three 
actuator  strips  are  patterned  in  series  with  the  specimen. 
Before  their  release  from  the  substrate,  both  actuator  and 
specimen  experience  the  same  tensile  strain  (residual  strain 
after  deposition).  However,  removal  of  the  substrate  also 
removes  the  uniform  strain  level  and  force  distribution  in 
the  freestanding  actuator  and  specimen,  which  is  re-organized 
according  to  their  respective  cross-sectional  areas.  Because  of 
its  larger  cross-sectional  area,  the  three  actuator  strips  (with 
the  length  La  and  cross-sectional  area  Aa )  will  apply  a  tensile 
force  on  the  specimen  (with  the  length  Ls  and  cross-sectional 
area  As).  The  amount  of  this  force  and  associated  strain  can  be 
controlled  by  the  actuator  design  utilizing  the  force  equilibrium 
condition 

F  =  ka8a  =  ks8s ,  (1) 

where  k  is  stiffness,  8  is  the  net  displacement  and  the  subscripts 
a  and  v  denote  the  actuator  and  specimen,  respectively.  The 
equilibrium  equation  can  thus  be  written  as  a  function  of  the 
length  (L),  cross-sectional  area  (A)  and  Young’s  modulus  (E) 
of  the  actuator  and  specimen 


where  n  is  the  number  of  actuator  strips.  Larger  value  of  n 
implies  higher  value  of  stress  and  strain  on  the  specimen. 
However,  experimental  results  deviate  more  than  25%  from 
the  model  for  n  >  5  because  of  the  compliance  of  the  outer 
actuator  strips  compared  to  the  inner  ones.  We  note  that 
the  total  displacement  in  the  fabricated  structure,  8 ,  is  the 
sum  of  the  actuator  and  specimen  displacement  (8a  +  8S  = 
8).  This  displacement  can  be  measured  by  fracturing  the 
specimen  after  the  experiment  and  then  measuring  the  gap 
between  the  fractured  surfaces.  Therefore,  irrespective  of  the 
overestimation  (if  at  all)  of  the  model,  the  resolution  of  strain 
measurement  is  governed  by  the  scanning  electron  microscope 
magnification.  Using  the  force  equilibrium  equation,  we  can 
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Silicon  nitride  deposition 


(a) 


Silicon  nitride  patterning  by 
plasma  etching 


Aluminum  (heater  line)  deposition 
by  lift-off  process  v. 


(d) 


Isotropic  etch  of  Si  substrate  using 
XeF2  to  make  beams  freestanding 


Figure  2.  Fabrication  processing  for  the  device  integrating  the  actuator,  specimen  and  heater  for  thermal  conductivity  characterization. 


relate  the  specimen  displacement  with  the  total  displacement 
as  follows: 


1  + 


EsAsLa 

nEaAaLs 


=  5. 


The  strain  in  the  specimen,  ss,  is  then  given  by 


&s_  _  SnEgAg 

LS  (fiEgAgLg  +  ESASLg ) 


(3) 

(4) 


For  this  study,  the  actuator  and  the  specimen  are  fabricated 
from  the  same  material.  The  specimen  strain  is  therefore  given 
by 


8nAa 

ss  = - .  (5) 

(nAaLs+AsLa) 

The  device  fabrication  process  is  schematically  described  in 
figure  2.  The  process  starts  with  a  deposition  of  50  nm  thick 
freestanding  LPCVD  Si3N4  using  ammonia  and  dichlorosilane 
at  820  °C  (figure  2(a)).  The  residual  stress  in  the  deposited 
film  is  about  1.2  GPa.  The  device  design  shown  in  figure  1(a) 
is  then  patterned  on  the  nitride  layer  using  lithography,  after 
which  the  exposed  nitride  is  etched  with  reactive  ion  etching 
using  CF4/O2  plasma.  This  is  shown  in  figure  2(b).  Using 
the  front  side  alignment,  negative  lithography  is  performed 
to  pattern  a  line  perpendicular  to  the  silicon  nitride  line  and 
175  nm  thick  aluminum  film  is  evaporated  to  complete  the 
lift-off  process.  Figure  1(c)  shows  the  aluminum  and  nitride 
lines  patterned  on  the  silicon  wafer,  which  are  still  attached  to 
the  substrate.  The  next  step  is  to  release  the  two  lines  using 
isotropic  XeF2  etching.  Figure  2(d)  shows  the  freestanding 
nitride  specimen  and  the  aluminum  heater  line  schematically, 
for  which  a  scanning  electron  micrograph  is  shown  in 
figure  1(b). 


3.  Experimental  procedure  and  analysis 

To  measure  thermal  conductivity  of  the  freestanding  nitride 
specimens,  we  pass  dc  current  through  the  aluminum  heaters. 
Since  the  heater  line  is  also  freestanding,  the  Joule  heating 
develops  a  parabolic  temperature  profile  with  the  highest 
temperature  in  the  middle  section,  which  intersects  the  silicon 
nitride  specimen  and  supplies  a  heat  flux  to  it.  In  a  typical 
experiment,  4-10  milliamps  of  current  is  passed  to  raise  the 
middle  section  of  the  heater  (and  the  silicon  nitride)  from  310 
to  375  K.  The  temperature  distribution  along  the  length  of  the 


Table  1.  Multi-physics  modules  and  input  parameters. 


Application 

mode 

Properties 

Silicon 

nitride 

Aluminum 

Conductive 

Resistivity  (12  m) 

1.72el2 

3.75e-8 

media  dc  mode 

Temperature 

0.0038 

0.003 

coefficient  of 
resistance  (1/K) 
Reference  temp  (K) 

300 

300 

Heat  transfer 

Thermal 

Varied  to 

237 

mode 

conductivity 

match 

(W  (m  K)-1) 

Density  (kg  m-3) 

measured 

temperature 

profile 

3100 

2700 

Heat  capacity, 

700 

904 

Solid  mechanics 

Cp  (J  (kg  K)-1) 
Young’s 

250e9 

70e9 

mode 

modulus  (Pa) 
Poisson’s  ratio 

0.23 

0.35 

Thermal  expansion 

2.3e-6 

23.1e-6 

coefficient  (1/K) 
Density  (kg  m-3) 

3100 

2700 

heater  is  measured  using  an  infrared  microscope  (Quantum 
Focus  Instruments,  MWIR-512  InSb  IR  FPA  camera)  with 
0.1  K  temperature  and  2  gm  spatial  resolution.  The  inset 
in  figure  1(b)  shows  a  thermal  micrograph  of  the  heater 
during  operation.  At  the  same  time,  the  temperature  profile 
along  the  length  of  the  silicon  nitride  specimen  is  measured. 
Calculation  of  the  thermal  conductivity  requires  the  heat 
flux  in  the  specimen,  which  can  be  approximated  from  the 
power  supply  to  the  heater,  or  more  accurately  by  performing 
finite  element  simulation  [35].  We  used  commercially 
available  software  COMSOL  Multiphysics™  to  carry  out  the 
finite  element  solution  of  the  described  problem.  Among  the 
various  application  modes,  electro-thermal  application  mode  is 
used  for  the  modeling  as  this  problem  incorporates  Joule 
heating  of  the  material.  This  mode  presents  the  coupling  of 
three  different  basic  application  modes  which  are  conductive 
media  dc  application,  heat  transfer  application  and  solid 
mechanics  application  modes.  Heat  loss  through  convection 
and  conduction  (through  the  air  to  the  etched  substrate  below 
the  specimen)  was  included  in  the  model.  Tables  1  and  2 
describe  primary  input  parameters  for  the  model.  The  physical 
properties  [37]  are  taken  for  bulk  scale  since  silicon  nitride  is 
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Boundary  2 


Specimen  -»  ? 


Boundary  4 


5xl0'5 


Boundary  3 


cq  y 


350 

340 

330 

320 

310 

300 


Figure  3.  Multi-physics  simulation  model  for  the  device  showing 
temperature  profile  at  zero  strain.  Infrared  microscopy  shows  the 
large  drop  in  thermal  conductivity  at  2.4%  strain  (inset). 


Table  2.  Geometry  parameters  for  the  aluminum  and  silicon  nitride 
structures. 


Material 

Length  (/xm) 

Width  (/xm) 

Thickness  (/x m) 

Aluminum 

35 

3 

0.175 

Silicon  nitride 

100 

5 

0.050 

an  amorphous  material  and  aluminum  film  is  too  thick  to  show 
any  size  effect. 

Figure  3  shows  the  finite  element  model,  where  the  four 
boundaries  are  numbered  for  the  ease  of  identification.  The 
conductive  media  dc  mode  solves  the  following  equation: 

_V.(aW  -Je)  =  Qj,  (6) 

where  a  is  the  electrical  conductivity,  V  is  the  potential,  JQ  is 
the  external  current  density  and  Qj  is  the  current  source.  This 
mode  assumes  that  the  model  has  symmetry  and  the  electric 
potential  varies  only  in  the  X  and  Y  directions  and  is  constant 
in  the  Z  (film  thickness)  direction.  For  Joule  heating,  the 
electrical  conductivity  is  temperature  dependent  and  maintains 
the  following  relationship: 


p0(l+a(T-T0))' 

where  pQ  is  the  resistivity  at  reference  temperature  T0  and  a 
is  the  temperature  coefficient  of  resistance.  One  end  of  the 
aluminum  heater  (boundary  3)  is  provided  with  an  electric 
potential  (V  =  V0)  and  the  opposite  end  (boundary  4) 
is  electrically  grounded,  as  shown  in  figure  4.  All  other 
boundaries  of  aluminum  and  silicon  nitride  are  chosen  to  be 
electrically  insulated  ( n.J  =  0). 

The  heat  transfer  mode  solves  the  following  governing 
equation: 

X(-kXT)  =  Q  +  qsT,  (8) 

where  k  is  the  thermal  conductivity,  XT  is  the  temperature 
gradient,  Q  is  the  heat  generation  and  qs  is  the 
production/absorption  coefficient.  The  four  ends  of  the 
aluminum  and  silicon  nitride  films  (boundaries  1,  2,  3  and  4) 
are  kept  at  constant  temperature  (T  =  To).  This  is  because  the 
thermal  mass  of  the  substrate  is  enormous  compared  to  the 
aluminum  and  nitride  films.  For  all  other  boundaries,  heat  flux 


boundary  condition  is  applied  through  the  following  equation: 

n.q  =  qQ  +  h(T0  -  T),  (9) 

where  n  is  the  outward  normal,  q  is  the  conductive  heat 
flux  vector,  q0  is  the  inward  heat  flux  from  external  sources 
(such  as  radiation)  and  h  is  the  convection  heat  transfer 
coefficient.  In  this  study  T0  is  taken  to  be  300  K.  Because  the 
experiments  were  conducted  at  lower  temperatures,  radiation 
was  neglected.  At  the  same  time,  heat  loss  through  air  (in 
this  case,  from  the  bottom  surface  of  the  nitride  specimen  to 
the  device  floor)  is  known  to  be  dominant  in  the  literature 
[38],  which  was  incorporated  through  an  ‘effective  convective 
heat  transfer  coefficient’ .  A  sensitivity  analysis  of  the  multi¬ 
physics  model  yielded  the  value  of  =  500  W  m-2  K-1  for 
the  etch  depth  of  about  15  gm  in  our  study,  for  which  the 
temperature  profiles  in  the  specimen  and  the  heater  were 
in  very  close  agreement  for  all  input  current  values.  A 
simple  scaling  analysis  may  suggest  higher  values  (h^  = 
1666  W  m-2  K-1),  but  is  unable  to  reproduce  the  actual 
temperature  profile  in  the  specimen. 

The  solid  mechanics  mode  employs  the  governing 
equation 

a=Es,  (10) 

where  a  is  the  stress,  E  is  Young’s  modulus  and  s  is  the  net 
strain.  The  net  strain  includes  applied  elastic  strain,  thermal 
strain  and  initial  strain.  Boundaries  1, 2, 3  and  4  are  constrained 
with  fixed  boundary  condition  (displacement  is  zero  in  all 
directions),  while  all  other  boundaries  are  considered  as 
constraint  free.  No  external  load  is  applied  to  the  boundaries. 

To  validate  the  model,  we  compare  the  simulated 
temperature  profile  along  the  aluminum  heater  line  with  that 
obtained  from  simulation  for  different  values  of  input  current. 
Figure  4  shows  the  excellent  agreement  between  the  actual  and 
simulated  heater  temperature  profiles.  It  also  shows  the  heat 
flux  value  from  the  simulation.  Since  the  governing  equations 
of  the  finite  element  model  do  not  have  any  strain-thermal 
conductivity  coupling,  we  first  measure  the  temperature  profile 
along  the  specimen  length  for  various  levels  of  mechanical 
strain.  To  extract  the  thermal  conductivity  of  the  specimen 
for  each  strain  values,  we  simulate  the  heat  flow  for  the 
calculated  heat  flux  values  for  various  trial  values  of  thermal 
conductivity  until  the  simulated  temperature  profile  closely 
matched  that  obtained  from  thermal  microscopy.  Figure  5 
shows  this  procedure  for  two  different  values  of  strain. 

4.  Results  and  discussion 

Figure  6  shows  the  experimental  results  for  room  temperature 
thermal  conductivity  of  50  nm  thick  freestanding  amorphous 
silicon  nitride  as  a  function  of  mechanical  strain  up  to  2.4%. 
The  fracture  strain  of  silicon  nitride  can  be  as  high  as  3% 
[39, 40] .  Also  shown  (in  latter  form)  are  the  data  taken  from  the 
literature  for  comparable  specimen  thicknesses.  The  literature 
data  are  shown  at  about  0.3%  strain  because  even  though 
the  specimens  are  freestanding,  their  clamped  edge  boundary 
conditions  do  not  relieve  the  residual  strain.  These  studies 
involve  SiNx  (also  known  as  low  stress  or  hydrogen  rich  silicon 
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Figure  4.  Comparison  of  the  experimental  and  simulation  results  on  the  heater  temperature  profile  as  a  function  of  input  current. 


Distance  from  heater,  microns 

Figure  5.  Comparison  of  experimental  and  simulation  results  on  the  temperature  profile  of  the  silicon  nitride  specimens  with  fitted  values  of 
thermal  conductivity. 


nitride)  with  x  varying  from  1  to  1.1,  for  which  the  residual 
strain  reported  is  about  0.2%-0.3%  [41].  When  interpolated 
between  0%  and  0.6%  strain,  our  results  are  in  excellent 
agreement  with  the  literature  except  for  two  instances  (studies 
D  and  E).  These  two  studies  involved  specimens  on  substrate, 
for  which  thermal  conductivity  is  expected  to  show  a  pseudo¬ 
size  effect  unless  corrected  for  thermal  contact  resistance.  This 
is  because  such  cross-plane  measurements  involve  nanoscale 
length  of  heat  flow  path  (L  =  specimen  thickness)  and  two 
interfacial  thermal  resistances  (pth)  [42] 

1  1  A  Pth 

- =  - +  — ,  (11) 

k(T,L)  k(T,L^oq)  L 


where  A  is  the  cross-section  of  the  specimen.  Lee  and  Cahill 
[33]  explained  their  observed  size  effect  with  the  above 
argument,  while  Bai  et  al  [36]  corrected  their  data  after 
developing  a  technique  to  measure  the  interfacial  thermal 
resistance. 

Two  concerns  related  to  the  experimental  technique 
deserve  further  scrutiny.  The  first  one  is  on  the  role  of 
interfacial  thermal  resistance,  which  is  typically  not  critical  for 
in-plane  thermal  conductivity  measurements.  This  is  verified 
by  plugging  in  representative  values  from  our  experiments 
in  equation  (11).  For  the  heater-specimen  interfacial  area 
A  =  25  x  1(T12  m2,  pth=  2-3  x  1(T8  m2K  W_1[33,  36] 
and  L  is  about  15  x  10-6  m),  the  size  dependent  term  is 
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Figure  6.  Thermal  conductivity  of  50  nm  thick  Si3N4  films  as  a  function  of  mechanical  strain.  The  letters  corresponding  to  0.3%  strain 
represent  data  from  the  literature.  The  symbol  (_L)  denotes  cross-plane  thermal  conductivity  data. 


only  about  6.25  x  10-14,  which  clearly  shows  that  interfacial 
thermal  resistance  is  not  a  critical  concern  for  this  study.  The 
second  concern  arises  from  the  amount  of  strain  developed 
in  the  aluminum  heater  due  to  the  strain  in  the  silicon 
nitride  specimen,  which  may  change  the  electrical  and  thermal 
properties  of  the  heater  and  introduce  experimental  errors.  For 
a  125  nm  thick  aluminum  film  on  silicon  oxide,  such  strain 
effects  have  been  shown  to  be  pronounced  [27],  which  calls  for 
a  detailed  analysis.  Since  the  freestanding  aluminum  heater  is 
attached  to  the  nitride  specimen  at  the  midpoint,  the  maximum 
strain  it  experiences  is  the  same  as  the  net  strain  in  the  width 
direction  of  the  nitride  specimen.  This  is  given  by 


&AI  —  £y 


( ~vcrx  +  ay) 
E 


^residual? 


(12) 


where  x  and  y  directions  are  along  the  length  and  width  of  the 
silicon  nitride  specimen,  respectively,  and  Poisson’s  ratio  v  = 
0.22.  The  residual  strain  is  about  0.4%  for  our  specimens  and 
is  bi-axial  and  tensile  in  nature.  The  actuation  method  applies 
uni-axial  tension  in  the  nitride  specimens  with  stress  ( ox )  about 
6  GPa  for  about  2.4%  strain  (the  fracture  stress  of  silicon  nitride 
varies  from  8  to  10  GPa).  For  these  values,  the  net  strain  in 
the  aluminum  heater  considering  the  Poisson  effect  is  —0.5%. 
In  other  words,  there  is  no  appreciable  strain  developed  in 
the  aluminum  heater  for  the  actuated  nitride  specimens.  This 
is  also  verified  by  our  measurements  of  electrical  resistance 
of  the  aluminum  heater  before  and  after  the  release  from 
substrate,  which  did  not  show  any  appreciable  change  in 
the  resistance  due  to  the  loading  of  the  heater  line.  For  zero 
strain-thermal  conductivity  measurements,  the  tensile  strain  in 
the  heater  segment  overlapping  the  nitride  specimen  is  about 
0.4%,  for  which  the  thermal  conductivity  was  assumed  to  be 
100  W(m  K)-1  [27]. 

To  explain  the  observed  influence  of  strain  on  thermal 
conductivity,  we  note  that  commonly  cited  size  effects  for  low¬ 
dimensional  materials  [23,  43]  and  hetero  structures  [44,  45] 
cannot  be  used  because  of  the  amorphous  microstructure  of  the 
specimens.  The  phonon  mean  free  path  of  amorphous  materials 
is  of  the  same  order  as  the  structural  disorder  [46],  rendering 


thermal  conductivity  size  independent.  Here,  the  phases  of  the 
atomic  vibrations  are  incoherent  and  the  oscillators  are  damped 
strongly  enough  to  pass  their  energy  to  the  neighbors  within 
half  the  period  of  oscillation;  yet  the  states  are  not  localized. 
In  other  words,  heat  transfer  in  amorphous  solids  is  a  random 
walk  of  energy  between  localized  oscillators  of  varying  sizes 
and  frequencies.  The  Cahill  and  Pohl  [47]  model  for  thermal 
conductivity  of  amorphous  materials  integrates  the  specific 
heat  (C),  phonon  velocity  ( v )  and  the  mean  free  path  (/)  over 
the  entire  range  of  frequencies 

1  d C 

km=~  —v(co)l(co)dto.  (13) 

3  Jo  da> 

In  real  amorphous  materials,  the  inherent  aperiodicity  gives 
rise  to  strong  localization  of  high  frequency  vibrations  or 
‘frac tons’.  Orbach  and  co-workers  therefore  suggested  an 
alternative  mode  of  heat  transfer,  modeled  by  hopping  energy 
through  highly  localized  vibrational  modes  (fractons)  that 
are  coupled  to  phonons  via  the  anharmonic  interaction  [48]. 
While  fractons  alone  in  a  harmonic  system  cannot  transfer 
heat,  it  has  been  shown  in  the  literature  that  anharmonicity 
induces  interaction  between  phonons  and  fractons  that  allow 
transfer  of  energy  from  a  localized  excitation  to  other.  Such 
fraction-phonon  interaction  resembles  spatial  overlapping  of 
the  localized  vibrational  modes  [49].  Without  such  overlap, 
amorphous  materials  would  exhibit  zero  thermal  conductivity. 
At  the  same  time,  the  overlap  must  be  minimal  because 
amorphous  materials  are  known  to  exhibit  very  low  thermal 
conductivity.  Therefore,  we  propose  that  overlap  is  extremely 
marginal  and  sensitive  to  tensile  deformation.  This  is  supported 
by  our  experimental  observation  of  about  one  order  of 
magnitude  reduction  in  thermal  conductivity  at  about  2.4% 
strain.  On  the  other  hand,  compressive  deformation  would 
only  increase  such  overlap,  which  will  eliminate  the  influence 
of  strain,  unless  the  stress  is  so  large  that  the  elastic  constants 
are  changed.  Interestingly,  this  conjecture  is  supported  by  both 
theory  and  recent  experiments  [26]. 

While  the  influence  of  strain  on  thermal  conductivity  is 
not  a  size  effect,  we  suggest  that  it  depends  on  the  vibrational 
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density  of  states,  which  in  turn  depends  on  atomic  bond  energy, 
distance  etc  parameters.  Unlike  crystalline  materials,  where  all 
the  modes  are  non-localized  (diffusive),  amorphous  materials 
show  the  spectra  of  both  localized  (fraction)  and  diffusive 
(phonon)  modes  [50].  Therefore,  the  thermal  conductivity  can 
be  expressed  as 

k(T)  =  km(T)  +  khov(T).  (14) 

Our  conceptual  model  suggests  that  the  strain-thermal 
conductivity  coupling  would  depend  on  the  relative 
contribution  of  these  two  modes,  which  in  turn,  depend  on 
the  lack  of  long-range  order  and  atomic  bond  energy. 
Crystalline  or  polycrystalline  materials  exemplify  diffusion- 
dominated  thermal  transport,  for  which  very  large  strain 
(~10%  [23])  or  stress  (~20  GPa  [51])  is  shown  to  produce 
only  modest  reduction  in  thermal  conductivity  at  the  extremely 
small  (<  10  nm)  length  scales.  The  lack  of  long-range  order  in 
amorphous  materials  gives  rise  to  the  predominance  of  fractons 
in  the  vibrational  density  of  states.  For  example,  the  sharpest 
peak  for  amorphous  silicon  is  at  the  frequency  of  66  meV, 
which  is  very  close  to  the  cutoff  (below  which  the  localization 
effect  is  too  strong  to  transport  heat)  value  of  72  meV  [52]. 
We  hypothesize  that  the  relative  contribution  of  fractons  in 
thermal  transport  of  amorphous  materials  is  dependent  on  the 
bond  energy  of  the  atoms.  While  the  vibrational  density  of 
states  of  silicon  nitride  is  not  available  in  the  literature,  we 
expect  the  peak  frequency  to  be  higher  than  amorphous  silicon 
because  the  Si-N  bond  energy  (78.32  kcal/mol)  is  much 
higher  than  the  Si-Si  bond  (53.31  kcal/mol)  [53].  Amorphous 
silicon  nitride  is  therefore  expected  to  show  stronger  vibration 
localization  compared  to  amorphous  silicon.  This  implies  that 
strain-thermal  conductivity  coupling  in  amorphous  solids  will 
also  depend  on  the  bond  type  and  energy.  Currently,  we  are 
performing  atomistic  simulations  to  verify  this  hypothesis. 

5.  Conclusion 

In  this  study,  we  investigated  the  concept  of  tuning  thermal 
conductivity  of  nanoscale  solids  using  externally  applied 
mechanical  strain.  We  developed  an  experimental  technique 
to  apply  controllable  strain  on  freestanding  amorphous  silicon 
nitride  specimens  and  measure  their  thermal  conductivity.  We 
demonstrate  the  technique  on  50  nm  thick  specimens  at  up  to 
2.4%  strain.  It  is  commonly  held  that  amorphous  solids  already 
show  minimum  values  of  thermal  conductivity  strain  at  zero 
strain  and  hence  are  expected  to  show  the  weakest  strain- 
transport  coupling.  In  contrast,  our  experimental  results  show 
very  strong  influence  on  thermal  conductivity,  which  decreased 
from  2.7  to  0.34  W  (m  K)-1  under  2.4%  strain.  We  propose  that 
thermal  conductivity  in  amorphous  solids  is  a  strong  function 
of  vibration  localization,  which  increases  with  tensile  strain. 
We  conclude  that  crystalline  solids  do  not  show  such  strong 
coupling  because  of  their  non-localized  vibrational  modes. 
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